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ABSTRACT 



Emission-lines in the form of filamentary structures is common in bright clusters characterized by short cooling times. In the Perseus 
cluster, cold molecular gas, tightly linked to the Ha filaments, has been recently revealed by CO observations. In order to understand 
the origin of these filamentary structures, we have investigated the evolution of the hot ICM gas perturbed by the AGN central 
activity in a Perseus like cluster. Using very-high resolution TreeSPH simulations combined with a multiphase model and a model of 
plasma bubbles, we have been able to follow the density and temperature evolution of the disturbed ICM gas around the bubbles. Our 
simulations show that a fraction of the 1-2 keV gas present at the center of clusters is trapped and entrained by the rising buoyant 
bubble to higher radius where the AGN heating is less efficient. The radiative cooling makes it cool in a few tens of Myr below 10'' K, 
forming cold filamentary structures in the wake and in the rim of the bubbles. The predicted cold gas formed outside the cluster center 
is in agreement with the total CO mass and density profile of the observed molecular gas as well as with the kinematics of the Ha 
filaments. This scenario explains the Ha and CO filaments observed in luminous clusters without contradicting the observed lack of 
1 keV gas. It also emphasizes that if the AGN feedback provides some heating (negative feedback) it also perturbs the ICM, increasing 
its cooling (positive feedback). 
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1. Introduction 

Cooling flows in clusters of galaxi es has been a puzzling 
topic for more than 30 years (Fabian| [l994l : iPeterson & Fabian 
2006). While early models of X-ray clusters jCowie & Binnev 
1977h applied to observations predict a large cold gas 
mass accretion rate, the recent data from Chandra and 
XMM-Newton ruled out such large amounts of cooling 
(iPeterson et alj l2003h . There are now numerous observa- 
tions showing that the AGN activity may be a clue to the 
cooling flow problem ('Fabian et al.' '2000*; 'McNamara et al 
2000; Owen et al. 2000; McNamar a et al 2000; Fabian et aT 
2003a b; Blanton et al. 20 031; IWise'et all (20061; IFabian et al 



. 12006- McNamara & NulsenFMOTi) . The 
that cooling flows 



are 



emerging picture is 
considerabl y reduced by the heat- 
ing provided by the AGN feedback (Briiggen & Kaiser '2001"; 
Churazov et al. 2001; OuiHs et al. 2001; Bohrinaeret al. 2002; 
Reynolds et al. 2002; Omm a et al.l200 4; Omma & Binnev 2004; 
Siiacki & SpringeL2006;, Heinz et al.ll2006;,Cattaneo & Tevssieri 



2007h . The i7ulk of the ICM gas remains hot (> 2keV) down to 
the cluster center, where it is accreted by the central black hole 
(Bondi accretion). 

However, in this picture, it is dif ficult to understand the ori- 
gin of t he observed cold ne i itral g as (lEdge'2001^, 'Edge & Fraveij 
I2003t ISalome & CombesI I2003L |2004; Salome et al. lOiM 
and the fila mentary Ha structures ( Conselice et alj I20()U 
ICrawford et al. 2005). The connection be tween radio source and 
cold CO gas ISalome & CombesI |2004 as well as the excess 
of optical/UV/IR light in clusters with central cooling time less 
than 3 Gyr are signs of the link between the filaments, the cool- 
ing and the AGN activity. Very deep Chandra observation of the 
Perseus cluster shows a clear correlation between the soft X-ray 
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emission (bet ween 0.5 and 1 keV) and the optical filamentary Ha 
structures (Fabian et alJl2003"bl l2006h . As the filaments cover a 
large range of cooling time within the hot gas, it is unlikely that 
they form by direct cooling of the 2 keV component. Moreover, 
the kinematics of the filament in Ha ((Hatch et al. 2006) and 
soft X-rays (Fabi an et al. [ (20061) matches the gas flow und er a 
buoyantly rising bubble. This led lFabian et"ai] ( l2003bl(2006 l) and 
iHatch et al] ( l2006l) to propose that the filaments have been pulled 
out from a central reservoir of cold and warm gas by the action 
of radio plasma bubbles inflated by AGN jets. 

In this letter, without answering the question of the inner 
quenching of the cooling flow, we present a series of very high 
resolution TreeSPH simulations of a Perseus type cluster includ- 
ing AGN feedback. We show that buoyantly rising bubbles en- 
train a fraction of the 1-2 ke V central ICM gas to much higher 
radius where it has time to cool and fall back to the cluster center, 
forming cold filamentary structures {T < 10"*). 



2. The model 

2.1. Initial mass distribution 

We have constructed a cluster model in order to fit the Perseus 
data from 'Sa nders et a n (12004V The initial conditions are based 
on two components, a dark halo and a hot gas phase. Both share 
the same j6-profile with /3 = 2/3 (pseudo-isothermal sphere), 
with a core radius of 40kpc and a central mass density of 
Pq - 1.38 X lO^Mokpc^l The gas represents 15% of the to- 
tal mass and is treated as self-gravitating in the fixed potential of 
the halo. Assuming a hydrogen mass fraction of 0.76, we found 
a central electron density of 0.073 e"/cm^ and a temperature of 
2.85 X 10'^ K. The model is truncated at a radius of 2 Mpc, giving 
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a total mass of 5.3 x 10'^ M©. With these parameters, the mean 
cooling time at the center is about 500 Myr. 

2.2. Gas physics 

Computing the interaction between the AGN and the hot ICM 
requires us to deal with a huge range of temperatures, spanning 
from more than 10^ K down to 10"* K. This may be achieved by 
assuming a multiphase behavior of the gas where hot and cold 
gas represents two distinct phases. 

The warm-hot gas phase (T > 10"* K) is modeled by an ideal 
gas with an adiabatic index of 5/3. The radiative cooling of the 
ga s is computed from the norm alized cooling function tabulated 
by [Sutherland & Dopitaj j 19931) . assuming an abundance of one 
third solar dSanders et al.l l2004l) . Hydrodynamical equations are 
integrated using the Lagrangian SPH technique. 

Below 10"* K, instead of quenching the cooling, the gas is 
assumed to change phase, becoming very cold and clumpy, 
as commonly observed in the ISM. In this phase, the semi- 
collisional nature of the gas is well modeled by the sticky parti- 
cle technique widely us ed to compute die behavior of the ISM 
in the disk of galaxies ('Schwarz"l98l'; 'Combes & Gerinll 19851; 
ISemelin & Co mbes 2002; Bournaud & Combes 200"^! 

We take into account the possibility of forming stars from 
the cold phase. This is computed following the standard pro- 
cess described in iKatz et al.l d 19961) which simulates quite well a 
Schmidt law. We also take into account supernova feedback by 
injecting energy in both thermal (10%) and kinetic (90%) form, 
assuming that 10'*^ erg are released for each solar mass formed. 



2.3. AGN feedback 

We have modeled the AGN feedback by instantaneously gen- 
erating bubbles of gas in the ICM This technique, which has 
already been used by s everal author s (Briiggen & Kaiseij 20011: 
Churazov et all 120011: 'Saxton et al.' 2001; iSiiacki & Springell 
2006; Siiacki et al. 112007; Gardini 2007), does not consider the 
evolution of the jet a t its origin that may be modeled by mass 
or 



where 



2001 



2003 



energ y injection ^Reynolds et al.' '2002; 'Ouilis et al 
Briiggen et alj 1200^ [B riiggen & Kaiser 2002; Briiggen 
Dalla Vecchia et al. 2004; Omma & Binney 2004; Omma et all 
2004; Heinz et al. 2006; Cattaneo & Teyssier lOQlh- We have 
implemented a new bubble generator allowing us to control 
precisely the relative pressure, density and specific energy of 
the bubble with respect to the ambient gas. We have also 
taken into account the presence of the puzzling cooler and 
denser gas f orming a brig ht shell around th e radio bubbles 
jFabian et al 1 2000; McNa mara et al.1 ,2000; Bl anton et al.l l200ll: 
LNulsen et al.' 2002; Blanton et al.' '2003; Fabian et all l2003bi 
.2006; McNamara & Nulsen 2007), asssuming that it is the result 
of ambient gas pushed out by the inflating bubble. From these 
assumptions, we can generate bubbles with densities differing 
from the ambient gas, without violating the total mass conser- 
vation. In a given sphere of radius Rb at a distance from 
the cluster center, the density p,, mass m,, specific energy m, and 
pressure P, of a particle / are modified as follows : 
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if R<Rs 
+ 1 ifR> Ru 



(2) 



(3) 



where the radius R is computed with respect to the bubble center 
The thickness of the rim is parametrized by R^, set to 1/4 of 
the bubble radius. The density and pressure of the bubble with 
respect of the ambiant cluster gas are controlled by parameters 
6o and q-q (see Table [TJ. When ao is set to 1, the pressure of 
the bubble is not modified with respect to the ICM (see Eq.[Tli. 
In this case, from first thermodynamic principles, no energy is 
injected into the system Q. The constant 6'q is computed in order 
to conserve the total mass. According to observations, the bubble 
rims are slightly denser and cooler than the ambient gas, while 
the pressure remains unchanged. 

2.4. Simulations and parameters 

The models have been run usi ng a modified v ersion of the par- 
allel TreeSPH code Gadget2 (ISpringelll200"5l) which conserves 
both energy and entropy. Additional physics have been included 
as described in paragraph l2.2l and l2.3l 

In order to reach a very high resolution in the cluster central 
region (up to 200 kpc), we have constructed our model using a 
multi-resolution technique. The models are made out of 7 cells 
where two consecutive cells differ by their resolution only by a 
factor of two, ensuring the stability of the integration of the mo- 
tion equations. This technique allows us to compute the evolu- 
tion of the bubbles within a central region of 200 kpc, reaching 
a mass resolution of 2.4 x 10* Mq and a spatial resolution of 
about 1 kpc, appreciably better than previous works using simi- 
lar Lagrangian codes. The total number of particles is 6'04r002, 
equivalent to 33'554'432 particles in a similar model where no 
multi-resolution technique is used. 

In order to explore the effect of the bubbles on the ambient 
gas, we have run several simulations exploring all the parame- 
ters. We have completed our study with a control run where no 
bubble is launched. This simulation is described in AppendixlAl 
Hereafter, we focus on a subset of 5 representative simula- 
tions. For each simulation, the parameters are listed in Table [T] 
Bubbles are created at T = 20 Myr and are followed for 600 Myr. 

Except for model 1, we have set ao - 1 (bubbles are in pres- 
sure equilibri um with the ambient gas). Th is choice is motivated 
by the work of lCattaneo & Tevssieii (|2P07') showing that the bulk 
of the AGN feedback energy can be deposited at the center of the 
cluster in the form of thermal energy. Consequently, as we have 
not implemented central thermal feedback, our models are sub- 
ject to the cooling catastrophe, occurring in the 10 central kpc. 



3. Results 

3. 1 . Global evolution of the bubbles 

The evolution of the bubbles for the 5 models is displayed in 
Fig. [1] The first stage is very similar to previous works. Once 



(1) 



' In fact the total energy is slightly modified, due to the particle 
mass and density changes which imply a change in the total potential. 
However, this change is very small compared to the total potential of 
the cluster and may be neglected. 
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Id 


Db [kpc] 


Rb [kpc] 




So 


Eb [10"" erg] 


rB[io'K] 


1 


50 


30 


1.4 


0.5 


0.8 


4.7 


2 


50 


30 


1.0 


0.5 


0.0 


4.1 


3 


50 


30 


1.0 


0.1 


0.0 


19.0 


4 


30 


10 


1.0 


0.5 


0.0 


4.5 


5 


10 


8 


1.0 


0.1 


0.0 


30.2 



Table 1. List of simulations and respective parameters. See 
Section|2]for the meaning of each parameter. 




Fig. 1. Temperature map of the evolution of bubbles for models 
1 (top) to 5 (bottom line). Colors represent the logarithm of the 
temperature between 10"* and 10*^ K as indicated by the color 
bar. Cold gas is represented by blue colors. The temperatures 
are computed from a thin slice of 2 kpc width centered on y = 0. 
The width and height of each box is 150 kpc. Time increases 
from left to right. 



generated, the bubbles quickly start moving upwards, reach- 
ing velocities of between 50 and 300 km/ s and forming typi- 
cal atomic mushrooms. However, contrary to Eulerian simula- 
tions, the latter phase does not show the formation of a torus. 
The bub ble is only weakly affected by hydrodynamic instabil- 
ities (see lAgertz et al.l (12006 ) for a complete discussion of the 
damping of Kelvin-Helmoltz and Rayleigh-Taylor instabilities 
in SPH codes) and dissolves only when reaching regions of sim- 
ilar specific entropy. 



3.2. Cold gas formation 

A common feature between the 5 simulations is the system- 
atic formation of very cold gas in the wake of the bubbles, far 
from the cluster cor^ In a few Myrs, a substantial fraction 
of gas cools from the ambient ICM temperature (1 - 2keV) 
down to 10'' K, with a small fraction of it reaching temper- 
atures below 10"* K and forming a thin trunk (see Fig. [TJ. 
The presence of cooler gas in the trunk of a buoyantly ris- 
ing bub ble already ha s been mentioned by different authors 
(Churazo vetaTI 120011: iBruggen & Kaiseil 1200 ll: ISaxton et all 
2001 :.Gardinill2007b . However, in al these simulations, the gas 
temperature never falls below 5 x 10^ K. The high resolution ob- 
tained in our simulations reinforced by the Lagrangian technique 
enables us to resolve higher densities, critical for the correct es- 
timation of the cooling time. 

In our simulations, the filaments are formed by 1 - 2keV 
gas (observed in the central regions of the Perseus cluster 
(iFabian et al.l l2006i) ) having short cooling time (between 400 
and 800 Myr) at the time the bubble is formed. In the case of 
large bubbles (models 1,2,3) this gas is situated below the bub- 
ble. When smaller bubbles are created near the center (models 
4,5), gas with sufficiently short cooling times is also found above 
the bubble. In both cases, this gas is trapped by the rising bubble 
and entrained outside the cooling cluster core where its suffi- 
ciently short cooling time allows it to cool quickly to very low 
temperatures. Fig. |2] shows the evolution of particles becoming 
part of the filaments (7 < 10'' K) at f = 400 Myr for model 2 and 
5 (See also Fig. IB. II and Fig. IB. 21 for the evolution of particles 
becoming part of the filaments at f = 600 Myr). At the beginning 
of simulation 2, particles being in the filaments at t - 400 Myr 
(resp. t - 600 Myr) are part of the bottom rim of the bubble 
at a distance of between 20 and 40 kpc (resp. 30 and 50 kpc) 
of the center, with a temperature between 1 and 2 ke V. During 
the rising of the bubble, this gas is dragged up to 80 kpc (resp. 
90 kpc). The radiative cooling slowly decreases the gas tempera- 
ture, which reaches 1 ke V around t - 300 Myr. At these temper- 
atures, the gas is no longer supported by the pressure and falls 
towards the center while being compressed. Its density is en- 
hanced by a factor of 10 which substantially decreases its cool- 
ing time 0. A very cold thin filament with a temperature below 
IO^'K is then formed. The main difference of model 5 is that 
the bubble is formed at smaller radius. Gas with a short enough 
cooling time (< 600 Myr) is initially situated above the bubble. 
As for the model 2, it is dragged at higher altitude (up to 60 kpc) 
and crossed by the bubble, before reaching low temperatures and 
forming the cold observed filament. 



3.3. Comparison witli observations 

While a more complete model, including thermal heat, visosity 
and turbulence will be needed in the future, in Appendix |C] we 
show that the CO gas observed in the Perseus cluster may be 
compatible with the production of cold gas by a dozen of bubble 
pairs similar to bubbles of model 5. The kinematics of the cool- 
ing gas (with temperature less than 10^ K along the filament is 
displayed in 

Fig. [3] The top of the filament is still entrained by the bubbles 
and has positive velocities between 100 and 200 km/ s, while the 
bottom of the filament is falling towards the center with negative 



- See Appendix IaI for an estimation of the cooling core in a control 
simulation run without feedback 

^ At these temperatures and densities, the cooling is so strong that it 
dominates the adiabatic compression heating. 
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Fig. 2. Evolution of the gas around the bubble of model 2 (top) 
and 5 (bottom) that forms the cold filament (r < 10'' K) at f = 
400 Myr. Colors have the same meaning as in Fig.[T] The plasma 
bubble is represented in grey. The box size is respectively 50x50 
(top) and 25 x 25 kpc(bottom). 
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Fig. 3. Velocity of gas in the filament of model 5 at f = 400 Myr 
mesured along the z axis, left : Distance of the filaments from 
the center, as a function of their speed, right : Surface density of 
the gas. The particles represented in the left part of the figure are 
contained in the red box. 

velocities. This result, showing that the filaments are stretching, 
is qual itatively in agreement with the observations o f Hatch et al.. 
(l2006h . 

4. Discussion and conclusion 

New very-high resolution SPH simulations of the Perseus cluster 
show that the cold gas observed in Ho' and CO filaments outside 
the cluster core is a natural product of the perturbations of the 
ICM by buoyantly rising bubbles. During the rise of a bubble, 
a fraction of the low cooling time ambient ICM gas around the 
bubble is entrained and cools at higher radius, forming cold gas 
below lO'^K in the trunk of the bubble. The st reaching of the 
filaments observed in Ha by dHatch et alj|2006l) is qualitatively 
reproduced. 



Because no central heating is taken into account, all our 
models are subject to the cooling catastrophe. Is the proposed 
scenario still valid in a model where the cooling flow is quenched 
by heating processes ? The role of the heating source is to 
suppress the formation of gas with a temperature lower than 
2 keV, but 1-2 k eV gas will still be present in abundance in 
the cluster center (iFabian et alJ 12006). As seen in section [3721 
this is precisely this gas that is trapped by the bubbles and later 
forms the cold filaments. The heating at the origin of the cool- 
ing flow quenching will thus have only a very weak influence 
on the amount of cold gas formed. In the new emerging picture 
of cooling flow, AGN feedback provides some heating (nega- 
tive feedback), but it is also responsible for a positive feedback, 
where ICM perturbations induce cold gas pr oduction. This is a 
step forward in th e 'cold feedback' scenario (' Soker et"aLll200ll : 
Pizzolato & Soker 2005h where the feedback takes place within 
a large region {R < 5 - 30kpc). 

Acknowledgements. We are grateful to the anonymous referee for interesting 
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Fig. A.l. Temperature map of the evolution of the cluster for the 
control simulation run with no bubbles. Colors represent the log- 
arithm of the temperature between lO"* and 10^ K as indicated 
by the color bar. Cold gas is then represented by blue colors. 
The temperatures are computed from a thin slice of 2 kpc width 
centered on y = 0. The width and height of each box is 150 kpc 




Fig. A.2. Cold gas surface density as a function of radius for the 
control run, at f = 400 (green) and t - 5 00 Myr (red). The das hed 
line corresponds to values derived from lSalome et al.l (l2006l) . 



Appendix A: Control run without feedback 

We present the result of the control simulation run without feed- 
back. Fig. lA.lK the equivalent of Fig. [T]i shows the temperature 
map of the cluster The cooling catastrophe occurs after 300 Myr. 
After 600 Myr, only the central 10 kpc have temperatures below 
10^ K. Fig. IA.3I shows the evolution of the mean temperature 
profile of the cluster The cluster core is defined as the region 
where the temperature drops below 10^ K. 

Fig. |A. 21 shows the cold gas surface density of the outlying 
regions {R > 10 kpc). Those values have been used to correct the 
equivalent values for models with bubble perturbation presented 
in in Appendix |C] 



Appendix B: Formation of cold filaments 

Fig. IB. 1 | and |B.2| show the evolution and interaction with the ris- 
ing bubble of the gas of model 2 and 5, that ends in the filament 
at f = 600 with a temperature below 10^ K. 




EO 40 



Fig. A.3. Evolution of the mean temperature profile of the ref- 
erence model, corresponding to a purely cooling flow (black 
: t = 20 Myr, blue : t = 400 Myr, red : f = 600 Myr). The 
dashed vertical line corresponds to our definition of the core ra- 
dius (T < 10^ K). Dots corresponds to the positions of particles 
forming a cold filament, at f = 400 and t = 600 Myr (black : 
/ - 20 Myr, blue : t - 400 Myr, red : t - 600 Myr, see also 
Fig.©. 
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Fig. B.l. Evolution of the gas around the bubble of model 2 that 
form the cold filament (T < lO^K) at f = 600 Myr. Colors 
represent the logarithm of the temperature between lO'' and 
10'^ K. The plasma bubble is represented in grey. The box size 
is 100x200 kpc. 
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Fig. B.2. Evolution of the gas around the bubble of model 4 that 
form the cold filament (r < lO^K) at f = 600 Myr. Colors 
represent the logarithm of the temperature between 10"* and 
10*^ K. The plasma bubble is represented in grey. The box size 
is 50 X 100 kpc. 



Appendix C: Comparison with CO observations 

In this appendix, we compute the amount of cold gas {T < 
10"* K) formed outside the cluster core {R > 10 kpc), due to 
the uplift of the ambient ICM gas during the bubble rise and 
compare it with available CO observation of the Perseus clus- 
ter. As no central heating is taken into account at the cluster 
center, a large amount of cold gas is naturally deposited within 
the central region (at a rate of about 16OM0/yr). From the 
amount of gas found for each simulation, we have then sub- 
tracted the equivalent value obtained from the control simulation 
(see Appendix lAb. where no bubbles have been launched. This 
ensures that the cold gas observed for R > 10 kpc is only due to 
ICM perturbation by the rising bubbles. 

Fig. IC.U compares the cold gas surface density of models 
1, 2, 3 and 5 0, at three different times, with the surface den- 
sity of the cold gas deduced from the Salome et al. (2006) data 
(dashed line). Clearly, all models fail to reproduce the amount 
of gas observed between 10 and 40 kpc. However, we have to 
recall that in those models, only two symmetric bubbles are re- 
sponsible for the ICM cooling outside the cluster core. In a more 
complete model, taking into account the recurrent or continuous 
activity of the AGN will generate more bubbles, multiplying the 
production of cold gas by a dex. It then means that models 1, 
2 and 3 probably over predict the cold gas, especially at a ra- 
dius larger than 40 kpc. However, the surface density profile of 
model 5 at f = 400 and f - 500 becomes remarkably similar to 
the Perseus observations. This argumentation is supported by the 
direct comparison of cold gas in the HorseShoe region (Pos 1 1 
in[Salome et al. (2007)). An amount of 0.9 ■ 10** Mq has been es- 
timated in this region, remarkably similar to the cold gas formed 
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Fig. C.l. Cold gas surface density as a function of radius for 
models 1, 2, 3, 5, at f = 300 (blue), t = 400 (green) and t = 
500 Myr (red). The dashed line corresponds to values derived 
from I Salome et al., (,2006,) . 



in the rim of the bubbles of model 2, 3 and 4, ranging between 
0.5 and lOx IO'^Mq. 

It is important to notice that a carefull comparison between 
the models and observations should also take into account the 
lifetime of the filaments with respect of the AGN duty cycle. 
Unfortunately, the lifetime of the filaments is difficult to pre- 
dict because it depends on several physical processes like heat 
conduction (Nipoti & Binnev 20(34i), viscosity or turbulence, not 
included in the present model. 



* As the cold gas production of model 4 is nearly two decades below 
the observations, we haven't displayed it 



